Tertiary-Quaternary subduction processes and related magmatism in the Alpine-Mediterranean region

During Tertiary to Quaternary times, convergence between Eurasia and Africa resulted in a variety of collisional orogens and different styles of subduction in the Alpine-Mediterranean region. Characteristic features of this area include arcuate orogenic belts and extensional basins, both of which can be explained by roll-back of subducted slabs and retreating subduction zones. After cessation of active subduction, slab detachment and post-collisional gravitational collapse of the overthickened lithosphere took place. This complex tectonic history was accompanied by the generation of a wide variety of magmas. Most of these magmas (e.g. low-K tholeiitic, calc-alkaline, shoshonitic and ultrapotassic types) have trace element and isotopic fingerprints that are commonly interpreted to reflect enrichment of their source regions by subduction-related fluids. Thus, they can be considered as ‘subduction-related’ magmas irrespective of their geodynamic relationships. Intraplate alkali basalts are also found in the region generally postdated the ‘subduction-related’ volcanism. These mantle-derived magmas have not been, or only slightly, influenced by subduction-related enrichment. This paper summarises the geodynamic setting of the Tertiary-Quaternary “subduction-related” magmatism in the different segments of the Alpine-Mediterranean region (Betic-Alboran-Rif province, Central Mediterranean, the Alps, Carpathian-Pannonian region, Dinarides and Hellenides, Aegean and Western Anatolia), and discusses the main characteristics and compositional variation of the magmatic rocks. Radiogenic and stable isotope data indicate the importance of continental crustal material in the genesis of these magmas. Interaction with crustal material probably occurred both in the upper mantle during subduction (‘source contamination’) and in the continental crust during ascent of mantle-derived magmas (either by mixing with crustal melts or by crustal contamination). The 87Sr/86Sr and 206Pb/204Pb isotope ratios indicate that an enriched mantle component, akin to the source of intraplate alkali mafic magmas along the Alpine foreland, played a key role in the petrogenesis of the ‘subduction-related’ magmas of the Alpine-Mediterranean region. This enriched mantle component could be related to mantle plumes or to long-term pollution (deflection of the central Atlantic plume and recycling of crustal material during subduction) of the shallow mantle beneath Europe since the late Mesozoic. In the first case, subduction processes could have had an influence in generating asthenospheric flow by deflecting nearby mantle plumes due to slab roll-back or slab break-off. In the second case, the variation in the chemical composition of the volcanic rocks in the Mediterranean region can be explained by “statistical sampling” of the strongly inhomogeneous mantle followed by variable degrees of crustal contamination

The interpretation of crustal dynamics data in terms of plate interactions and active tectonics of the Anatolian plate and surrounding regions in the Middle East

A detailed study was made of the consequences of the Arabian plate convergence against Eurasia and its effects on the tectonics of Anatolia and surrounding regions of the eastern Mediterranean. A primary source of information is time rates of change of baseline lengths and relative heights determined by repeated SLR measurements. These SLR observations are augmented by a network of GPS stations in Anatolia, Aegea, and Greece, established and twice surveyed since 1988. The existing SLR and GPS networks provide the spatial resolution necessary to reveal the details of ongoing tectonic processes in this area of continental collision. The effort has involved examining the state of stress in the lithosphere and relative plate motions as revealed by these space based geodetic measurements, seismicity, and earthquake mechanisms as well as the aseismic deformations of the plates from conventional geodetic data and geological evidence. These observations are used to constrain theoretical calculations of the relative effects of: (1) the push of the Arabian plate; (2) high topography of Eastern Anatolia; (3) the geometry and properties of African-Eurasian plate boundary; (4) subduction under the Hellenic Arc and southwestern Turkey; and (5) internal deformation and rotation of the Anatolian plate

Active tectonics of the Hellenic subduction zone

This thesis is remarkable for the wide range of the techniques and observations used and for its insights, which cross several disciplines. It begins by solving a famous puzzle of the ancient world, which is what was responsible for the tsunami that destroyed settlements in the eastern Mediterranean in 365 AD. By radiocarbon dating of preserved marine organisms, Shaw demonstrates that the whole of western Crete was lifted out of the sea by up to 10 metres in a massive earthquake at that time, which occured on a previously unknown fault. The author shows that the resulting tsunami would have the characteristics described by ancient writers, and uses modern GPS measurements and coastiline geomorphology to show that the strain build-up near Crete requires such a tsunami-earthquake about every 6.000 years - a major insight into Mediterranean tsunami hazard. A detailed seismological study of earthquakes in the Cretan arc over the last 50 years reveals other important features of its behaviour that were previously unknown. Finally, she provides fundamental insights into the limitations of radiocarbon dating marine organisms, relating to how they secrete carbon into their skeletons. The thesis resulted in three major papers in top journals

Global correlations between maximum magnitudes of subduction zone interface thrust earthquakes and physical parameters of subduction zones

Peer reviewedPublisher PD

Constraining the long-term evolution of the slip rate for a major extensional fault system in the central Aegean, Greece, using thermochronology

The brittle/ductile transition is a major rheologic boundary in the crust yet little is known about how or if rates of tectonic processes are influenced by this boundary. In this study we examine the slip history of the large-scale Naxos/Paros extensional fault system (NPEFS), Cyclades, Greece, by comparing published slip rates for the ductile crust with new thermochronological constraints on slip rates in the brittle regime. Based on apatite and zircon fission-track (AFT and ZFT) and (U–Th)/He dating we observe variable slip rates across the brittle/ductile transition on Naxos. ZFT and AFT ages range from 11.8 ± 0.8 to 9.7 ± 0.8 Ma and 11.2 ± 1.6 to 8.2 ± 1.2 Ma and (U–Th)/He zircon and apatite ages are between 10.4 ± 0.4 to 9.2 ± 0.3 Ma and 10.7 ± 1.0 to 8.9 ± 0.6 Ma, respectively. On Paros, ZFT and AFT ages range from 13.1 ± 1.4 Ma to 11.1 ± 1.0 Ma and 12.7 ± 2.8 Ma to 10.5 ± 2.0 Ma while the (U–Th)/He zircon ages are slightly younger between 8.3 ± 0.4 Ma and 9.8 ± 0.3 Ma. All ages consistently decrease northwards in the direction of hanging wall transport. Most of our new thermochronological results and associated thermal modeling more strongly support the scenario of an identical fault dip and a constant or slightly accelerating slip rate of 6–8 km Myr− 1 on the NPEFS across the brittle/ductile transition. Even the intrusion of a large granodiorite body into the narrowing fault zone at 12 Ma on Naxos does not seem to have affected the thermal structure of the area in a way that would significantly disturb the slip rate. The data also show that the NPEFS accomplished a minimum total offset of 50 km between 16 and 8 Ma

Slip distribution of the 2017 M(w)6.6 Bodrum-Kos earthquake: resolving the ambiguity of fault geometry

SUMMARY The 2017 July 20, Mw6.6 Bodrum–Kos earthquake occurred in the Gulf of Gökova in the SE Aegean, a region characterized by N–S extension in the backarc of the easternmost Hellenic Trench. The dip direction of the fault that ruptured during the earthquake has been a matter of controversy where both north- and south-dipping fault planes were used to model the coseismic slip in previous studies. Here, we use seismic (seismicity, main shock modelling, aftershock relocations and aftershock mechanisms using regional body and surface waves), geodetic (GPS, InSAR) and structural observations to estimate the location, and the dip direction of the fault that ruptured during the 2017 earthquake, and the relationship of this event to regional tectonics. We consider both dip directions and systematically search for the best-fitting locations for the north- and south-dipping fault planes. Comparing the best-fitting planes for both dip directions in terms of their misfit to the geodetic data, proximity to the hypocenter location and Coulomb stress changes at the aftershock locations, we conclude that the 2017 earthquake ruptured a north-dipping fault. We find that the earthquake occurred on a 20–25 km long, ∼E–W striking, 40° north-dipping, pure normal fault with slip primarily confined between 3 and 15 km depth, and the largest slip exceeding 2 m between depths of 4 and 10 km. The coseismic fault, not mapped previously, projects to the surface within the western Gulf, and partly serves both to widen the Gulf and separate Kos Island from the Bodrum Peninsula of SW Anatolia. The coseismic fault may be an extension of a mapped, north-dipping normal fault along the south side of the Gulf of Gökova. While all of the larger aftershocks are consistent with N–S extension, their spatially dispersed pattern attests to the high degree of crustal fracturing within the basin, due to rapid trenchward extension and anticlockwise rotation within the southeastern Aegean

Segmented Hellenic slab rollback driving Aegean deformation and seismicity

The NE dipping slab of the Hellenic subduction is imaged in unprecedented detail using teleseismic receiver function analysis on a dense 2-D seismic array. Mapping of slab geometry for over 300 km along strike and down to 100 km depth reveals a segmentation into dipping panels by along-dip faults. Resolved intermediate-depth seismicity commonly attributed to dehydration embrittlement is shown to be clustered along these faults. Large earthquakes occurrence within the upper and lower plate and at the interplate megathrust boundary show a striking correlation with the slab faults suggesting high mechanical coupling between the two plates. Our results imply that the general slab rollback occurs here in a differential piecewise manner imposing its specific stress and deformation pattern onto the overriding Aegean plate

Flexural deformation of the continental lithosphere

Prior work focused primarily on the Adriatic and northern Ionian regions. The results of these studies have been summarized previously, and so are only briefly discussed. More recent work focuses on two different topics: (1) analysis of foredeep basin geometry, sedimentary style, and thrust belt structure in light of the kinematics at the associated plate boundary and subduction zone dynamics; and (2) the evolution and plate strength of early Proterozoic lithosphere